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We present a search for standard model (SM) Higgs boson to WW^ production in dilepton plus missing 
transverse energy final states in data collected by the CDF II detector corresponding to 4.8 fb — 1 of integrated 
luminosity. To maximize sensitivity, the multivariate discriminants used to separate signal from background in 
the opposite-sign dilepton event sample are independently optimized for final states with zero, one, or two or 
more identified jets. All significant Higgs boson production modes (gluon fusion, associated production with 
either a W or Z boson, and vector boson fusion) are considered in determining potential signal contributions. We 
also incorporate a separate analysis of the same-sign dilepton event sample which potentially contains additional 
signal events originating from associated Higgs boson production mechanisms. Cross section limits relative to 
the combined SM predictions are presented for a range of Higgs boson mass hypotheses between 110 and 200 
GeV/c 2 . 
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1. Introduction 

The Higgs boson is introduced into the SM to ex- 
plain electroweak symmetry breaking and the origins 
of particle mass. This paper presents an analysis 
searching for a Higgs boson decaying to VFW^*) , which 
is the dominant decay mode for mjj > 135 GeV/c 2 
The small cross section of the dominant gluon fusion 
production mechanism ((JNNLiAgg — » H) = 0.439 pb 
for rriH — 160 GeV/c 2 @, d, Q) makes observation 
of the signal difficult within the hadron collider envi- 
ronment. This analysis searches for H — ■> WW^-*' — » 
£ + £~vv, where I = e, /i, or r with final states e + e~, 
e^/i^, anc j ^+ y- _ j n addition, a search is made for 
associated Higgs boson production in events with two 
like-sign dileptons, with final states e ± e ± , e ± /x ± , and 
H ± ti ± . The data corresponds to 4.8 fb _1 of integrated 
luminosity collected at the Fermilab Tevatron in pp 
collisions at y/s = 1.96 TeV. 



2. Event Selection 

The CDF II detector is described in detail in 
Ref. [f|. The geometry of the detector is character- 
ized by the azimuthal angle <p and the pseudorapidity 
rj = — ln[tan((9/2)], where 9 is the polar angle mea- 
sured from the proton beam direction. The transverse 
energy Et = EsinO, where E is the energy in the 
calorimeter towers associated with a cluster of energy 
deposition. Transverse momentum pt is the track 
momentum component transverse to the beam-line. 
The missing transverse energy vector $ T is defined 
as — ^\ E T n T , where n l T is the transverse compo- 
nent of the unit vector pointing from the interaction 
point to the energy deposition in calorimeter tower 
i. This is corrected for the px of muons, which do 
not deposit all of their energy in the calorimeter, and 
tracks which point to uninstrumented regions in the 
calorimeter. The missing transverse energy ]fi T is de- 
fined as \$ T \. Strongly interacting partons produced 



in the pp collision undergo fragmentation that results 
in highly collimated jets of hadronic particles. Jet can- 
didates are reconstructed using the calorimeter towers 
with corrections to improve the estimated energy Q 
and are required to have Et > 15 GeV and \r/\ < 2.5. 

The tlvv candidates are selected from two opposite- 
sign (OS) leptons. At least one lepton is required 
to satisfy the real-time trigger selection and have 
Pt > 20 GeV/c. This requirement is loosened to 10 
GeV/c for the second lepton to increase Higgs boson 
kinematic acceptance. The z-positions of the leptons 
in a candidate at the point of closest approach to the 
beam-line are required to be within 4 cm of each other 
to reduce backgrounds from overlapping pp collisions. 

There are several sources of background with the 
same or a similar final state: WW — > tlvv, ZZ — ► 
tlvv or lUl, WZ -> lllv or IvvV, it -> bbllvv, Drell- 
Yan (DY) Z — > It where the measured large lfi T is due 
to resolution tails, where the photon converts to 
an e~e + pair, and VF+jets where a jet is misiden- 
tificd as a lepton. To reduce the contribution of W7 
and VF+jets, all leptons must be isolated with require- 
ments on both the surrounding Et in the calorimeters 
and pt of surrounding tracks. Candidates are further 
required to have a dilepton invariant mass Mu > 16 
GeV/c 2 to suppress IF7, and exactly 2 leptons to sup- 
press contributions with additional leptons. 

To reduce the DY background, we require 
#T spC c > 25 GeV (reduced to # Tspcc > 15 GeV for 
electron-muon events where DY background is inher- 
ently smaller), where 



J Tspcc 



Ip T sin( A(f> 



> - 

T> nearest 2 

< - 

nearest 2 



The Achjf is the azimuthal separation between 

SpT , nearest 

the Iprp and the momentum vector of the nearest lep- 
ton or jet candidate. With this definition, we require 
the missing energy transverse to the nearest lepton 
or jet in the event to be greater than the minimum 
threshold if $ T points along that object. Thus a mis- 
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measurement of the energy of one lepton or jet candi- 
date would not allow the event to enter the sample. 

We consider individually final states with zero, one, 
or two or more jet candidates in the event. This al- 
lows us to tune our multivariate discriminants on the 
different relative sizes of signal and background con- 
tributions within each jet multiplicity bin. The dom- 
inant background in the zero jet bin is WW, while 
in the one jet bin the background from DY is of a 
similar size as the WW background. The dominant 
background in the two or more jets bin is it produc- 
tion. To suppress the it contribution, all events with 
two or more jets containing one or more jets with a 
tight, secondary vertex 6-tag are rejected. After all se- 
lection requirements we observe (summing over all jet 
multiplicities) a total of 1467 candidate events com- 
pared against an expectation of 1479±141 background 
events and 25 ± 3 signal events for a SM Higgs boson 
with a mass of 165 GeV/c 2 . 

2.1. Event Selection for Same Sign 
Analysis 

Like-sign, or same-sign (SS), dileptons occur in 
VH — > VWW production when the vector boson (Z 
or W) and one of the W bosons from the Higgs bo- 
son decays leptonically. The dominant backgrounds to 
this search are from the charge mismeasurement of a 
real lepton, photon conversion to e~e + , and misiden- 
tification of a jet as a lepton. 

The event selection for OS events is modified to 
select two like-sign leptons and reduce potential back- 
grounds. The SS channel does not use forward elec- 
tron candidates, which have a high rate of charge mis- 
measurement. To reduce the numbers of conversion 
electrons and jets misidentified as leptons, the pr re- 
quirement for the second lepton is increased from 10 
to 20 GeV/c. As the decay of the third boson most 
often results in the production of additional jets, we 
also require one or more jet candidates in the final 
state. Like-sign VH events tend to have lower values 
of Ip T because the neutrinos are not necessarily back- 
to-back. To increase the signal acceptance, no $>T apcc 
requirement is made. After all selection requirements 
we observe a total of 64 candidate events compared 
against an expectation of 62 ± 11 background events 
and 1.7 ± 0.2 signal events for a SM Higgs boson with 
a mass of 165 GeV/c 2 . 



3. Data Modeling 

The geometric and kinematic acceptance for the 
WW, WZ, ZZ, W-y, DY, ti, and all of the signal 
processes (gg — > H, WH , ZH , VBF) are determined 
using a Monte Carlo calculation of the collision fol- 
lowed by a GEANT3-based simulation of the CDF II 



detector response 0j- The Monte Carlo generator 
used for WW is MC@NLO 3, while for WZ, ZZ, 
DY, it, and the signal processes pythia [9-] is used, 
and Wj is modeled with the generator described in 
Ref. fldl. We use the CTEQ5L [llj parton distribution 
functions (PDF) to model the momentum distribution 
of the initial-state partons. 

The expected yields of background processes are 
normalized to cross sections calculated at partial 
next-to next-to-leading order (ttjl2|), next-to-leading 
order (WW, WZ, and ZZ [l3(), or leading-order 
with estimated higher-order corrections (W7 [l0| and 
DY 14]). The gluon fusion cross section has been 
calculated to next-to next-to-leading logarithmic ac- 
curacy @, H, H[ . Associated production cross sections 
at next-to next-to-leading order (NNLO) are obtained 
from Ref. [l5j]. VBF cross sections at next-to- leading 
order (NLO) are obtained from Ref. [l6j . 

The background from VF+jets is estimated from a 
sample of data events with one identified lepton and 
one jet that fulfills loose isolation requirements and 
contains a track or energy cluster similar to those re- 
quired for lepton identification. The contribution of 
each event to the total yield is scaled by the probabil- 
ity that such a jet is identified as a lepton. This prob- 
ability is determined from multijet events collected in 
independent jet-triggered data samples. A correction 
is applied for the small real lepton contribution using 
single W and Z boson Monte Carlo simulation. 

4. Cross Checks 

Several control regions are constructed to validate 
different aspects of the data modeling. The sample 
sizes of the control regions are designed to be large 
enough to give statistically meaningful tests. These 
control regions, described below, show good agree- 
ment between the background model predictions and 
the observed distributions in data. 

One of the most fundamental control regions is the 
Drell-Yan sample selected using the default OS can- 
didate selection in Sec. [5] with a reversed $ T < 25 
GeV (15 GeV in the case of epi events) requirement 
and the additional requirement that the dilepton in- 
variant mass fall in the Z mass region. This control 
sample provides a high statistics test of acceptance, 
lepton identification efficiency, and trigger efficiency 
calculations. 

A W+jets control sample is constructed from events 
which satisfy the SS candidate selection in Sec. l2.1l but 
have zero jet candidates. This control sample contains 
contributions primarily from jets misidentified as lep- 
tons, with additional contributions from conversion 
electrons and charge mismeasurement of a real lepton, 
and checks the modeling of misidentified particles. 

A it control sample is constructed using the default 
OS candidate selection in Sec. [3 with a requirement 
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of two or more jets in the event, and the 6-tag veto 
requirement reversed such that only events with one 
or more jets with a 6-tag are accepted. 

A second DY control sample tests the $ T model, 
and is constructed using the default OS candidate se- 
lection in Sec. [2] with the exception of requiring 15 
< #T sp cc < 25 GeV instead of # Tspoc > 25 GeV (in 
the case of e/j, events, 10 < $t SP cc < 1<> GeV instead 
of -$Tspcc > 15 GeV). No requirements are made on 
the jet multiplicity or the dilepton invariant mass. 



leptons' pt, the $ T , the $ T significance ($ T / \/TJE~t) , 
the fij- and the A(j>jp 

spec T ^//p,nearest 

Neural network outputs are shown in Figure [1] for 
each channel at a Higgs boson mass of 165 GeV/c 2 . 
Plots of the input variables and neural network out- 
puts for all Higgs boson masses are available in 
Ref. [3. 



6. Systematics 



5. Analysis Technique 

For each search channel, a NeuroBayes® [TtJ neu- 
ral network is trained on a weighted combination of 
signal and background events from Monte Carlo inde- 
pendently for each Higgs boson mass hypothesis (14 
total). Each neural network has three layers consisting 
of input nodes, hidden nodes, and one output node. 
Once the neural network is trained, templates are cre- 
ated for all signal and background processes. For the 
OS zero and one jet channels, "high S/B" and "low 
S/B" templates based on the signal-to-background ra- 
tio of the different dilepton combinations are consid- 
ered separately. The templates are used as the final 
discriminant in calculating the 95% credibility level 
(C.L.) limits. 

The OS zero jet channel uses 5 input variables. The 
inputs are: the AR between the leptons, the A<fi be- 
tween the leptons, the scalar sum of the leptons' px 
and the $ T , the likelihood ratio for Higgs to WW pro- 
duction (LRhww), and the likelihood ratio for WW 
production (LRww)- The likelihood ratios are de- 
termined using a matrix element method which cal- 
culates the probability for each event to have been 
produced by several relevant SM processes. 

The OS one jet channel uses 8 input variables. The 
inputs are: the Mu, the transverse mass of the sum 
of the leptons' 4- momenta and ($ T , MetX, MetY, 0) 
where MetX and MetY are the x and y components of 
the Jfjj-i the AR between the leptons, the scalar sum 
of the leptons' pr and jets' Et and the $ T (Ht), the 
-^Tspcc' Ihe leading lepton px, the subleading lepton 
Pt, and the leading lepton energy. 

The OS two or more jets channel uses 8 input vari- 
ables. The inputs are: the Mu, the leading lepton pt, 
the subleading lepton pt, the Ht, the AR between the 
leptons, the Acf> between the leptons, the A<p between 
the vector sum of the leptons' pr and the $J T , and the 
vector sum of the first and second jet Et- 

The SS one or more jets channel uses 13 input vari- 
ables. The inputs are: the Mu, the leading lepton en- 
ergy, the leading lepton pt, the subleading lepton pt, 
the Ht, the leading jet Et, the A</> between the lep- 
tons, the vector sum of all jets' Et, the vector sum of 



Systematic uncertainties associated with the Monte 
Carlo simulation affect the signal, WW, WZ, ZZ, 
Wj, DY, and it acceptances taken from the simu- 
lated event samples. Uncertainties originating from 
lepton selection and trigger efficiency measurements 
are propagated through the acceptance calculation, 
giving uncertainties typically around 2%, but as high 
as 7% for the different signal and background pro- 
cesses. For the SS channel, there is an additional sys- 
tematic for the simulated backgrounds which enter the 
sample only when a lepton's charge is mismeasured 
(WW, tt, and DY) which is taken as half the differ- 
ence between the mismeasurement rates found in DY 
events and those found in Monte Carlo. 

We also assign an acceptance uncertainty due to po- 
tential contributions from higher-order effects. In the 
case of WW we take half of the difference between the 
leading-order (PYTHlA-based [§|) and next-to-leading 
order (mc@NLO [|j]) acceptances. The other processes 
are only simulated at leading-order and for these 
modes we assign the full difference observed in WW, 
leading to a 10% uncertainty. The largest uncertainty 
on the Drell-Yan process originates from modeling of 
the fake $ T in these events. Based on discrepancies 
in the DY Jft T control region, we assign a 20% uncer- 
tainty in the zero jet channel and 25% uncertainties in 
the one and two or more jets channels. For processes 
that produce no final state jets at leading-order (WW, 
Wj, and DY), we also assign a jet modeling system- 
atic to account for uncertainites in the Monte Carlo 
modeling of the pt boost and additional jet produc- 
tion from initial and final state radiation. The uncer- 
tainties for these effects are anti-correlated between 
jet channels since the effect of extra jet production 
is to move events from one channel to another. For 
Monte Carlo samples not currently generated over the 
entire data run range, we take additional uncertain- 
ties on the acceptance for the corresponding signal or 
background process. This uncertainty comes from the 
observed difference in the leading-order WW accep- 
tance for the corresponding partial run range versus 
that for the full run range. The acceptance variations 
due to PDF model uncertainties is assessed to be on 
the order of 2% using the 20 pairs of PDF sets de- 
scribed in Ref. [l9l |. 

For the Wj background contribution, there is an 
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Figure 1: Neural network outputs for each channel (OS zero and one jet for both high and low S/B, OS two or more 
jets, and SS one or more jets) at mjj = 165 GeV/c 2 . The Higgs boson signal is shown at 10 times the SM prediction. 



additional uncertainty of 20% from the detector ma- 
terial description and conversion veto efficiency. The 
systematic uncertainty on the VF+jets background 
contribution is determined from differences in the 
measured probability that a jet is identified as a lepton 
for jets collected using different jet Er trigger thresh- 
olds. These variations correspond to changing the par- 
ton composition of the jets and the relative amount of 
contamination from real leptons. 

The cross section calculation uncertainties are 6% 
on diboson production, 10% on it and W-f produc- 
tion, and 5% on Drell-Yan production. All signal and 
background estimates obtained from simulation have 



an additional 6% uncertainty originating from the lu- 
minosity measurement [20] . 

Most systematic uncertainties on the signal pro- 
cesses are assessed using the same techniques de- 
scribed for the background processes. Uncertainties 
on the theoretical cross sections vary by Higgs boson 
production mechanism. Associated production cross 
sections are known to NNLO, so the theoretical un- 
certainty on these cross sections is small, less than 
5% [l6|. VBF production is known only to NLO, 
so the residual theoretical uncertainty is higher, on 
the order of 10% (lq |. Gluon fusion is a QCD pro- 
cess; thus the theoretical uncertainty is still signifi- 
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cant even though it is known to NNLO. We use the 
hnnlo program [2l|, [22J, which computes the theo- 
retical cross section for gg — ► H production at the 
Tevatron to NNLO, to assign theoretical uncertain- 
ties due to changing the renormalization and factor- 
ization scales (11%) and the gluon PDF model (5%). 
We again use the hnnlo program to evaluate accep- 
tance uncertainties for gg — > H production for both 
changes in scale and the PDF model. To evaluate 
the change in acceptance, we reweight the Higgs bo- 
son pt and 77 distributions obtained from PYTHIA to 
match the HNNLO calculations. Reweighting events 
based on Higgs boson pt changes the relative accep- 
tances for the jet channels, while reweighting events 
based on Higgs boson rj primarily affects the lepton 
acceptance. 



7. Results 

To determine the sensitivity we construct a binned 
likelihood using the six neural network templates 
shown in Figure [TJ two templates each for the OS 
zero and one jet channels (high S/B and low S/B), 
one template for the OS two or more jets channel, and 
one template for the SS one or more jets channel. All 
the signal and background normalizations are allowed 
to float, but ratios of signal and background contribu- 
tions are constrained to their expectations with a set 
of Gaussian constraints determined from the system- 
atic uncertainties. The total signal yield is allowed to 
float. 

The 95% C.L. limits are determined with a set 
of 50,000 Monte Carlo background-only experiments 
based on expected yields varied within the assigned 
systematics. For each experiment a test statistic is 
formed from the difference in the likelihood value for 
the background-only model versus that for the signal 
plus background model. The median expected 95% 
C.L. limit at a Higgs mass of 165 GeV/c 2 is 1.3±qI 
times the SM prediction, while the observed limit is 
1.3 times the SM prediction. Results for 14 values of 
run are shown in Table Q] and Figure [H 
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Figure 2: The ratio of expected and observed 95% C.L. production cross section limits to SM predictions for 14 different 
Higgs boson masses for the combination of all opposite-sign and same-sign H — > WW'*' decay channels. Error bands 
are shown on the expected limit, while the observed limit is shown as a solid line. 



Table I The ratio of expected and observed 95% C.L. production cross section limits to SM predictions for 14 different 
Higgs boson masses for the combination of all opposite-sign and same-sign H — > W 7 !^'*' decay channels. 
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